Ancient and recent origins of shared polymorphisms in yeast.

Shared genetic polymorphisms between populations and species can be ascribed to ancestral variation or to more recent gene flow. Here, we mapped shared polymorphisms in Saccharomyces cerevisiae and its sister species Saccharomyces paradoxus, which diverged 4-6 million years ago. We used a dense map of single-nucleotide diagnostic markers (mean distance 15.6 base pairs) in 1,673 sequenced S. cerevisiae isolates to catalogue 3,852 sequence blocks (≥5 consecutive markers) introgressed from S. paradoxus, with most being recent and clade-specific. The highly diverged wild Chinese S. cerevisiae lineages were depleted of introgressed blocks but retained an excess of individual ancestral polymorphisms derived from incomplete lineage sorting, perhaps due to less dramatic population bottlenecks. In the non-Chinese S. cerevisiae lineages, we inferred major hybridization events and detected cases of overlapping introgressed blocks across distinct clades due to either shared histories or convergent evolution. We experimentally engineered, in otherwise isogenic backgrounds, the introgressed PAD1-FDC1 gene pair that independently arose in two S. cerevisiae clades and revealed that it increases resistance against diverse antifungal drugs. Overall, our study retraces the histories of divergence and secondary contacts across S. cerevisiae and S. paradoxus populations and unveils a functional outcome.

Genome instability footprint under rapamycin and hydroxyurea treatments.

The mutational processes dictating the accumulation of mutations in genomes are shaped by genetic background, environment and their interactions. Accurate quantification of mutation rates and spectra under drugs has important implications in disease treatment. Here, we used whole-genome sequencing and time-resolved growth phenotyping of yeast mutation accumulation lines to give a detailed view of the mutagenic effects of rapamycin and hydroxyurea on the genome and cell growth. Mutation rates depended on the genetic backgrounds but were only marginally affected by rapamycin. As a remarkable exception, rapamycin treatment was associated with frequent chromosome XII amplifications, which compensated for rapamycin induced rDNA repeat contraction on this chromosome and served to maintain rDNA content homeostasis and fitness. In hydroxyurea, a wide range of mutation rates were elevated regardless of the genetic backgrounds, with a particularly high occurrence of aneuploidy that associated with dramatic fitness loss. Hydroxyurea also induced a high T-to-G and low C-to-A transversion rate that reversed the common G/C-to-A/T bias in yeast and gave rise to a broad range of structural variants, including mtDNA deletions. The hydroxyurea mutation footprint was consistent with the activation of error-prone DNA polymerase activities and non-homologues end joining repair pathways. Taken together, our study provides an in-depth view of mutation rates and signatures in rapamycin and hydroxyurea and their impact on cell fitness, which brings insights for assessing their chronic effects on genome integrity.

Chemical-genomic profiling identifies genes that protect yeast from aluminium, gallium, and indium toxicity.

Aluminium, gallium, and indium are group 13 metals with similar chemical and physical properties. While aluminium is one of the most abundant elements in the Earth's crust, gallium and indium are present only in trace amounts. However, the increased use of the latter metals in novel technologies may result in increased human and environmental exposure. There is mounting evidence that these metals are toxic, but the underlying mechanisms remain poorly understood. Likewise, little is known about how cells protect themselves from these metals. Aluminium, gallium, and indium are relatively insoluble at neutral pH, and here we show that they precipitate in yeast culture medium at acidic pH as metal-phosphate species. Despite this, the dissolved metal concentrations are sufficient to induce toxicity in the yeast Saccharomyces cerevisiae. By chemical-genomic profiling of the S. cerevisiae gene deletion collection, we identified genes that maintain growth in the presence of the three metals. We found both shared and metal-specific genes that confer resistance. The shared gene products included functions related to calcium metabolism and Ire1/Hac1-mediated protection. Metal-specific gene products included functions in vesicle-mediated transport and autophagy for aluminium, protein folding and phospholipid metabolism for gallium, and chorismate metabolic processes for indium. Many of the identified yeast genes have human orthologues involved in disease processes. Thus, similar protective mechanisms may act in yeast and humans. The protective functions identified in this study provide a basis for further investigations into toxicity and resistance mechanisms in yeast, plants, and humans.

Predictive evolution of metabolic phenotypes using model-designed environments.

Adaptive evolution under controlled laboratory conditions has been highly effective in selecting organisms with beneficial phenotypes such as stress tolerance. The evolution route is particularly attractive when the organisms are either difficult to engineer or the genetic basis of the phenotype is complex. However, many desired traits, like metabolite secretion, have been inaccessible to adaptive selection due to their trade-off with cell growth. Here, we utilize genome-scale metabolic models to design nutrient environments for selecting lineages with enhanced metabolite secretion. To overcome the growth-secretion trade-off, we identify environments wherein growth becomes correlated with a secondary trait termed tacking trait. The latter is selected to be coupled with the desired trait in the application environment where the trait manifestation is required. Thus, adaptive evolution in the model-designed selection environment and subsequent return to the application environment is predicted to enhance the desired trait. We experimentally validate this strategy by evolving Saccharomyces cerevisiae for increased secretion of aroma compounds, and confirm the predicted flux-rerouting using genomic, transcriptomic, and proteomic analyses. Overall, model-designed selection environments open new opportunities for predictive evolution.

Adaptation of the yeast gene knockout collection is near-perfectly predicted by fitness and diminishing return epistasis.

Adaptive evolution of clonally dividing cells and microbes is the ultimate cause of cancer and infectious diseases. The possibility of constraining the adaptation of cell populations, by inhibiting proteins enhancing the evolvability, has therefore attracted interest. However, our current understanding of how genes influence adaptation kinetics is limited, partly because accurately measuring adaptation for many cell populations is challenging. We used a high-throughput adaptive laboratory evolution platform to track the adaptation of >18,000 cell populations corresponding to single-gene deletion strains in the haploid yeast deletion collection. We report that the preadaptation fitness of gene knockouts near-perfectly (R2= 0.91) predicts their adaptation to arsenic, leaving at the most a marginal role for dedicated evolvability gene functions. We tracked the adaptation of another >23,000 gene knockout populations to a diverse range of selection pressures and generalized the almost perfect (R2=0.72-0.98) capacity of preadaptation fitness to predict adaptation. We also reconstructed mutations in FPS1, ASK10, and ARR3, which together account for almost all arsenic adaptation in wild-type cells, in gene deletions covering a broad fitness range and show that the predictability of arsenic adaptation can be understood as a by global epistasis, where excluding arsenic is more beneficial to arsenic unfit cells. The paucity of genes with a meaningful evolvability effect on adaptation diminishes the prospects of developing adjuvant drugs aiming to slow antimicrobial and chemotherapy resistance.

Genetically controlled mtDNA deletions prevent ROS damage by arresting oxidative phosphorylation.

Deletion of mitochondrial DNA in eukaryotes is currently attributed to rare accidental events associated with mitochondrial replication or repair of double-strand breaks. We report the discovery that yeast cells arrest harmful intramitochondrial superoxide production by shutting down respiration through genetically controlled deletion of mitochondrial oxidative phosphorylation genes. We show that this process critically involves the antioxidant enzyme superoxide dismutase 2 and two-way mitochondrial-nuclear communication through Rtg2 and Rtg3. While mitochondrial DNA homeostasis is rapidly restored after cessation of a short-term superoxide stress, long-term stress causes maladaptive persistence of the deletion process, leading to complete annihilation of the cellular pool of intact mitochondrial genomes and irrevocable loss of respiratory ability. This shows that oxidative stress-induced mitochondrial impairment may be under strict regulatory control. If the results extend to human cells, the results may prove to be of etiological as well as therapeutic importance with regard to age-related mitochondrial impairment and disease.

Genome-Wide Association Study Reveals Host Factors Affecting Conjugation in Escherichia coli.

The emergence and dissemination of antibiotic resistance threaten the treatment of common bacterial infections. Resistance genes are often encoded on conjugative elements, which can be horizontally transferred to diverse bacteria. In order to delay conjugative transfer of resistance genes, more information is needed on the genetic determinants promoting conjugation. Here, we focus on which bacterial host factors in the donor assist transfer of conjugative plasmids. We introduced the broad-host-range plasmid pKJK10 into a diverse collection of 113 Escherichia coli strains and measured by flow cytometry how effectively each strain transfers its plasmid to a fixed E. coli recipient. Differences in conjugation efficiency of up to 2.7 and 3.8 orders of magnitude were observed after mating for 24 h and 48 h, respectively. These differences were linked to the underlying donor strain genetic variants in genome-wide association studies, thereby identifying candidate genes involved in conjugation. We confirmed the role of fliF, fliK, kefB and ucpA in the donor ability of conjugative elements by validating defects in the conjugation efficiency of the corresponding lab strain single-gene deletion mutants. Based on the known cellular functions of these genes, we suggest that the motility and the energy supply, the intracellular pH or salinity of the donor affect the efficiency of plasmid transfer. Overall, this work advances the search for targets for the development of conjugation inhibitors, which can be administered alongside antibiotics to more effectively treat bacterial infections.

Domestication reprogrammed the budding yeast life cycle.

Domestication of plants and animals is the foundation for feeding the world human population but can profoundly alter the biology of the domesticated species. Here we investigated the effect of domestication on one of our prime model organisms, the yeast Saccharomyces cerevisiae, at a species-wide level. We tracked the capacity for sexual and asexual reproduction and the chronological life span across a global collection of 1,011 genome-sequenced yeast isolates and found a remarkable dichotomy between domesticated and wild strains. Domestication had systematically enhanced fermentative and reduced respiratory asexual growth, altered the tolerance to many stresses and abolished or impaired the sexual life cycle. The chronological life span remained largely unaffected by domestication and was instead dictated by clade-specific evolution. We traced the genetic origins of the yeast domestication syndrome using genome-wide association analysis and genetic engineering and disclosed causative effects of aneuploidy, gene presence/absence variations, copy number variations and single-nucleotide polymorphisms. Overall, we propose domestication to be the most dramatic event in budding yeast evolution, raising questions about how much domestication has distorted our understanding of the natural biology of this key model species.

Aborting meiosis allows recombination in sterile diploid yeast hybrids.

Hybrids between diverged lineages contain novel genetic combinations but an impaired meiosis often makes them evolutionary dead ends. Here, we explore to what extent an aborted meiosis followed by a return-to-growth (RTG) promotes recombination across a panel of 20 Saccharomyces cerevisiae and S. paradoxus diploid hybrids with different genomic structures and levels of sterility. Genome analyses of 275 clones reveal that RTG promotes recombination and generates extensive regions of loss-of-heterozygosity in sterile hybrids with either a defective meiosis or a heavily rearranged karyotype, whereas RTG recombination is reduced by high sequence divergence between parental subgenomes. The RTG recombination preferentially arises in regions with low local heterozygosity and near meiotic recombination hotspots. The loss-of-heterozygosity has a profound impact on sexual and asexual fitness, and enables genetic mapping of phenotypic differences in sterile lineages where linkage analysis would fail. We propose that RTG gives sterile yeast hybrids access to a natural route for genome recombination and adaptation.

Machine Learning Prediction of Resistance to Subinhibitory Antimicrobial Concentrations from Escherichia coli Genomes.

Escherichia coli is an important cause of bacterial infections worldwide, with multidrug-resistant strains incurring substantial costs on human lives. Besides therapeutic concentrations of antimicrobials in health care settings, the presence of subinhibitory antimicrobial residues in the environment and in clinics selects for antimicrobial resistance (AMR), but the underlying genetic repertoire is less well understood. Here, we used machine learning to predict the population doubling time and cell growth yield of 1,407 genetically diverse E. coli strains expanding under exposure to three subinhibitory concentrations of six classes of antimicrobials from single-nucleotide genetic variants, accessory gene variation, and the presence of known AMR genes. We predicted cell growth yields in the held-out test data with an average correlation (Spearman's ρ) of 0.63 (0.36 to 0.81 across concentrations) and cell doubling times with an average correlation of 0.59 (0.32 to 0.92 across concentrations), with moderate increases in sample size unlikely to improve predictions further. This finding points to the remaining missing heritability of growth under antimicrobial exposure being explained by effects that are too rare or weak to be captured unless sample size is dramatically increased, or by effects other than those conferred by the presence of individual single-nucleotide polymorphisms (SNPs) and genes. Predictions based on whole-genome information were generally superior to those based only on known AMR genes and were accurate for AMR resistance at therapeutic concentrations. We pinpointed genes and SNPs determining the predicted growth and thereby recapitulated many known AMR determinants. Finally, we estimated the effect sizes of resistance genes across the entire collection of strains, disclosing the growth effects for known resistance genes in each individual strain. Our results underscore the potential of predictive modeling of growth patterns from genomic data under subinhibitory concentrations of antimicrobials, although the remaining missing heritability poses a challenge for achieving the accuracy and precision required for clinical use. IMPORTANCE Predicting bacterial growth from genome sequences is important for a rapid characterization of strains in clinical diagnostics and to disclose candidate novel targets for anti-infective drugs. Previous studies have dissected the relationship between bacterial growth and genotype in mutant libraries for laboratory strains, yet no study so far has examined the predictive power of genome sequence in natural strains. In this study, we used a high-throughput phenotypic assay to measure the growth of a systematic collection of natural Escherichia coli strains and then employed machine learning models to predict bacterial growth from genomic data under nontherapeutic subinhibitory concentrations of antimicrobials that are common in nonclinical settings. We found a moderate to strong correlation between predicted and actual values for the different collected data sets. Moreover, we observed that the known resistance genes are still effective at sublethal concentrations, pointing to clinical implications of these concentrations.

Mutagenic mechanisms of cancer-associated DNA polymerase ϵ alleles.

A single amino acid residue change in the exonuclease domain of human DNA polymerase ϵ, P286R, is associated with the development of colorectal cancers, and has been shown to impart a mutator phenotype. The corresponding Pol ϵ allele in the yeast Saccharomyces cerevisiae (pol2-P301R), was found to drive greater mutagenesis than an entirely exonuclease-deficient Pol ϵ (pol2-4), an unexpected phenotype of ultra-mutagenesis. By studying the impact on mutation frequency, type, replication-strand bias, and sequence context, we show that ultra-mutagenesis is commonly observed in yeast cells carrying a range of cancer-associated Pol ϵ exonuclease domain alleles. Similarities between mutations generated by these alleles and those generated in pol2-4 cells indicate a shared mechanism of mutagenesis that yields a mutation pattern similar to cancer Signature 14. Comparison of POL2 ultra-mutator with pol2-M644G, a mutant in the polymerase domain decreasing Pol ϵ fidelity, revealed unexpected analogies in the sequence context and strand bias of mutations. Analysis of mutational patterns unique to exonuclease domain mutant cells suggests that backtracking of the polymerase, when the mismatched primer end cannot be accommodated in the proofreading domain, results in the observed insertions and T>A mutations in specific sequence contexts.

Genomic Epidemiology and Evolution of Escherichia coli in Wild Animals in Mexico.

Escherichia coli is a common bacterial species in the gastrointestinal tracts of warm-blooded animals and humans. Pathogenicity and antimicrobial resistance in E. coli may emerge via host switching from animal reservoirs. Despite its potential clinical importance, knowledge of the population structure of commensal E. coli within wild hosts and the epidemiological links between E. coli in nonhuman hosts and E. coli in humans is still scarce. In this study, we analyzed the whole-genome sequencing data of a collection of 119 commensal E. coli strains recovered from the guts of 55 mammal and bird species in Mexico and Venezuela in the 1990s. We observed low concordance between the population structures of E. coli isolates colonizing wild animals and the phylogeny, taxonomy, and ecological and physiological attributes of the host species, with distantly related E. coli strains often colonizing the same or similar host species and distantly related host species often hosting closely related E. coli strains. We found no evidence for recent transmission of E. coli genomes from wild animals to either domesticated animals or humans. However, multiple livestock- and human-related virulence factor genes were present in E. coli of wild animals, including virulence factors characteristic of Shiga toxin-producing E. coli (STEC) and atypical enteropathogenic E. coli (aEPEC), where several isolates from wild hosts harbored the locus of enterocyte effacement (LEE) pathogenicity island. Moreover, E. coli isolates from wild animal hosts often harbored known antibiotic resistance determinants, including those against ciprofloxacin, aminoglycosides, tetracyclines, and beta-lactams, with some determinants present in multiple, distantly related E. coli lineages colonizing very different host animals. We conclude that genome pools of E. coli colonizing the guts of wild animals and humans share virulence and antibiotic resistance genes, underscoring the idea that wild animals could serve as reservoirs for E. coli pathogenicity in human and livestock infections.IMPORTANCEEscherichia coli is a clinically important bacterial species implicated in human- and livestock-associated infections worldwide. The bacterium is known to reside in the guts of humans, livestock, and wild animals. Although wild animals are recognized as potential reservoirs for pathogenic E. coli strains, the knowledge of the population structure of E. coli in wild hosts is still scarce. In this study, we used fine resolution of whole-genome sequencing to provide novel insights into the evolution of E. coli genomes from a small yet diverse collection of strains recovered within a broad range of wild animal species (including mammals and birds), the coevolution of E. coli strains with their hosts, and the genetics of pathogenicity of E. coli strains in wild hosts in Mexico. Our results provide evidence for the clinical importance of wild animals as reservoirs for pathogenic strains and highlight the need to include nonhuman hosts in the surveillance programs for E. coli infections.

A High-Throughput Method for Screening for Genes Controlling Bacterial Conjugation of Antibiotic Resistance.

The rapid horizontal transmission of antibiotic resistance genes on conjugative plasmids between bacterial host cells is a major cause of the accelerating antibiotic resistance crisis. There are currently no experimental platforms for fast and cost-efficient screening of genetic effects on antibiotic resistance transmission by conjugation, which prevents understanding and targeting conjugation. We introduce a novel experimental framework to screen for conjugation-based horizontal transmission of antibiotic resistance between >60,000 pairs of cell populations in parallel. Plasmid-carrying donor strains are constructed in high-throughput. We then mix the resistance plasmid-carrying donors with recipients in a design where only transconjugants can reproduce, measure growth in dense intervals, and extract transmission times as the growth lag. As proof-of-principle, we exhaustively explore chromosomal genes controlling F-plasmid donation within Escherichia coli populations, by screening the Keio deletion collection in high replication. We recover all seven known chromosomal gene mutants affecting conjugation as donors and identify many novel mutants, all of which diminish antibiotic resistance transmission. We validate nine of the novel genes' effects in liquid mating assays and complement one of the novel genes' effect on conjugation (rseA). The new framework holds great potential for exhaustive disclosing of candidate targets for helper drugs that delay resistance development in patients and societies and improve the longevity of current and future antibiotics. Further, the platform can easily be adapted to explore interspecies conjugation, plasmid-borne factors, and experimental evolution and be used for rapid construction of strains.IMPORTANCE The rapid transmission of antibiotic resistance genes on conjugative plasmids between bacterial host cells is a major cause of the accelerating antibiotic resistance crisis. There are currently no experimental platforms for fast and cost-efficient screening of genetic effects on antibiotic resistance transmission by conjugation, which prevents understanding and targeting conjugation. We introduce a novel experimental framework to screen for conjugation-based horizontal transmission of antibiotic resistance between >60,000 pairs of cell populations in parallel. As proof-of-principle, we exhaustively explore chromosomal genes controlling F-plasmid donation within E. coli populations. We recover all previously known and many novel chromosomal gene mutants that affect conjugation efficiency. The new framework holds great potential for rapid screening of compounds that decrease transmission. Further, the platform can easily be adapted to explore interspecies conjugation, plasmid-borne factors, and experimental evolution and be used for rapid construction of strains.

A yeast living ancestor reveals the origin of genomic introgressions.

Genome introgressions drive evolution across the animal1, plant2 and fungal3 kingdoms. Introgressions initiate from archaic admixtures followed by repeated backcrossing to one parental species. However, how introgressions arise in reproductively isolated species, such as yeast4, has remained unclear. Here we identify a clonal descendant of the ancestral yeast hybrid that founded the extant Saccharomyces cerevisiae Alpechin lineage5, which carries abundant Saccharomyces paradoxus introgressions. We show that this clonal descendant, hereafter defined as a 'living ancestor', retained the ancestral genome structure of the first-generation hybrid with contiguous S. cerevisiae and S. paradoxus subgenomes. The ancestral first-generation hybrid underwent catastrophic genomic instability through more than a hundred mitotic recombination events, mainly manifesting as homozygous genome blocks generated by loss of heterozygosity. These homozygous sequence blocks rescue hybrid fertility by restoring meiotic recombination and are the direct origins of the introgressions present in the Alpechin lineage. We suggest a plausible route for introgression evolution through the reconstruction of extinct stages and propose that genome instability allows hybrids to overcome reproductive isolation and enables introgressions to emerge.

Protein kinase A controls yeast growth in visible light.

BACKGROUND: A wide variety of photosynthetic and non-photosynthetic species sense and respond to light, having developed protective mechanisms to adapt to damaging effects on DNA and proteins. While the biology of UV light-induced damage has been well studied, cellular responses to stress from visible light (400-700 nm) remain poorly understood despite being a regular part of the life cycle of many organisms. Here, we developed a high-throughput method for measuring growth under visible light stress and used it to screen for light sensitivity in the yeast gene deletion collection. RESULTS: We found genes involved in HOG pathway signaling, RNA polymerase II transcription, translation, diphthamide modifications of the translational elongation factor eEF2, and the oxidative stress response to be required for light resistance. Reduced nuclear localization of the transcription factor Msn2 and lower glycogen accumulation indicated higher protein kinase A (cAMP-dependent protein kinase, PKA) activity in many light-sensitive gene deletion strains. We therefore used an ectopic fluorescent PKA reporter and mutants with constitutively altered PKA activity to show that repression of PKA is essential for resistance to visible light. CONCLUSION: We conclude that yeast photobiology is multifaceted and that protein kinase A plays a key role in the ability of cells to grow upon visible light exposure. We propose that visible light impacts on the biology and evolution of many non-photosynthetic organisms and have practical implications for how organisms are studied in the laboratory, with or without illumination.

Intragenic repeat expansion in the cell wall protein gene HPF1 controls yeast chronological aging.

Aging varies among individuals due to both genetics and environment, but the underlying molecular mechanisms remain largely unknown. Using a highly recombined Saccharomyces cerevisiae population, we found 30 distinct quantitative trait loci (QTLs) that control chronological life span (CLS) in calorie-rich and calorie-restricted environments and under rapamycin exposure. Calorie restriction and rapamycin extended life span in virtually all genotypes but through different genetic variants. We tracked the two major QTLs to the cell wall glycoprotein genes FLO11 and HPF1 We found that massive expansion of intragenic tandem repeats within the N-terminal domain of HPF1 was sufficient to cause pronounced life span shortening. Life span impairment by HPF1 was buffered by rapamycin but not by calorie restriction. The HPF1 repeat expansion shifted yeast cells from a sedentary to a buoyant state, thereby increasing their exposure to surrounding oxygen. The higher oxygenation altered methionine, lipid, and purine metabolism, and inhibited quiescence, which explains the life span shortening. We conclude that fast-evolving intragenic repeat expansions can fundamentally change the relationship between cells and their environment with profound effects on cellular lifestyle and longevity.

Shared Molecular Targets Confer Resistance over Short and Long Evolutionary Timescales.

Pre-existing and de novo genetic variants can both drive adaptation to environmental changes, but their relative contributions and interplay remain poorly understood. Here we investigated the evolutionary dynamics in drug-treated yeast populations with different levels of pre-existing variation by experimental evolution coupled with time-resolved sequencing and phenotyping. We found a doubling of pre-existing variation alone boosts the adaptation by 64.1% and 51.5% in hydroxyurea and rapamycin, respectively. The causative pre-existing and de novo variants were selected on shared targets: RNR4 in hydroxyurea and TOR1, TOR2 in rapamycin. Interestingly, the pre-existing and de novo TOR variants map to different functional domains and act via distinct mechanisms. The pre-existing TOR variants from two domesticated strains exhibited opposite rapamycin resistance effects, reflecting lineage-specific functional divergence. This study provides a dynamic view on how pre-existing and de novo variants interactively drive adaptation and deepens our understanding of clonally evolving populations.

Inhibiting conjugation as a tool in the fight against antibiotic resistance.

Antibiotic resistance, especially in gram-negative bacteria, is spreading globally and rapidly. Development of new antibiotics lags behind; therefore, novel approaches to the problem of antibiotic resistance are sorely needed and this commentary highlights one relatively unexplored target for drug development: conjugation. Conjugation is a common mechanism of horizontal gene transfer in bacteria that is instrumental in the spread of antibiotic resistance among bacteria. Most resistance genes are found on mobile genetic elements and primarily spread by conjugation. Furthermore, conjugative elements can act as a reservoir to maintain antibiotic resistance in the bacterial population even in the absence of antibiotic selection. Thus, conjugation can spread antibiotic resistance quickly between bacteria of the microbiome and pathogens when selective pressure (antibiotics) is introduced. Potential drug targets include the plasmid-encoded conjugation system and the host-encoded proteins important for conjugation. Ideally, a conjugation inhibitor will be used alongside antibiotics to prevent the spread of resistance to or within pathogens while not acting as a growth inhibitor itself. Inhibiting conjugation will be an important addition to our arsenal of strategies to combat the antibiotic resistance crisis, allowing us to extend the usefulness of antibiotics.

Prediction of antibiotic resistance in Escherichia coli from large-scale pan-genome data.

The emergence of microbial antibiotic resistance is a global health threat. In clinical settings, the key to controlling spread of resistant strains is accurate and rapid detection. As traditional culture-based methods are time consuming, genetic approaches have recently been developed for this task. The detection of antibiotic resistance is typically made by measuring a few known determinants previously identified from genome sequencing, and thus requires the prior knowledge of its biological mechanisms. To overcome this limitation, we employed machine learning models to predict resistance to 11 compounds across four classes of antibiotics from existing and novel whole genome sequences of 1936 E. coli strains. We considered a range of methods, and examined population structure, isolation year, gene content, and polymorphism information as predictors. Gradient boosted decision trees consistently outperformed alternative models with an average accuracy of 0.91 on held-out data (range 0.81-0.97). While the best models most frequently employed gene content, an average accuracy score of 0.79 could be obtained using population structure information alone. Single nucleotide variation data were less useful, and significantly improved prediction only for two antibiotics, including ciprofloxacin. These results demonstrate that antibiotic resistance in E. coli can be accurately predicted from whole genome sequences without a priori knowledge of mechanisms, and that both genomic and epidemiological data can be informative. This paves way to integrating machine learning approaches into diagnostic tools in the clinic.

Clonal Heterogeneity Influences the Fate of New Adaptive Mutations.

The joint contribution of pre-existing and de novo genetic variation to clonal adaptation is poorly understood but essential to designing successful antimicrobial or cancer therapies. To address this, we evolve genetically diverse populations of budding yeast, S. cerevisiae, consisting of diploid cells with unique haplotype combinations. We study the asexual evolution of these populations under selective inhibition with chemotherapeutic drugs by time-resolved whole-genome sequencing and phenotyping. All populations undergo clonal expansions driven by de novo mutations but remain genetically and phenotypically diverse. The clones exhibit widespread genomic instability, rendering recessive de novo mutations homozygous and refining pre-existing variation. Finally, we decompose the fitness contributions of pre-existing and de novo mutations by creating a large recombinant library of adaptive mutations in an ensemble of genetic backgrounds. Both pre-existing and de novo mutations substantially contribute to fitness, and the relative fitness of pre-existing variants sets a selective threshold for new adaptive mutations.